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9 Suppose that Pn{x) is the nth Lagrange polynomial that agrees 
with the function f at the distinct numbers xq, xi , . . . , x^. 

• Although this polynomial is unique, there are alternate algebraic 
representations that are useful in certain situations. 

• The divided differences of f with respect to xq, xi , . . . , Xn are used 
to express Pn{x) in the form 

Pn(x) = ao+a^ (x-Xo)+a2(x-xo)(x-xi) + - • •+an(x-xo) • • • {x-Xn-^) 

for appropriate constants ao,a-\, . . . ,an. 
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written in the form of the above equation, then evaluating Pn(x) at 
xq leaves only the constant term ao; that is, 

ao = PniXo) = f{Xo) 

9 Similarly, when P(x) is evaluated at xi, the only nonzero terms in 
the evaluation of Pn(xi) are the constant and linear terms, 

f{Xo) + ai (Xi - Xo) = Pn(Xi ) = f(Xi ) 
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A Definition 



• The zeroth divided difference of the function f with respect to x/, 
denoted f[x/], is simply the value of f at x,: 

f[Xi] = m 

• The remaining divided differences are defined recursively. 
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• The first divided difference of f with respect to x, and x/+i is 
denoted f[x/,x/+i] and defined as 

X/+1 — X/ 

• The second divided difference, f[x,,x,+i,x,+2], is defined as 

.r^ ^ ^ 1 _ f[Xi+^,Xi+2]-f[Xi,Xi+^] 

i[Xi, X/+1 , X/+2J - 

X/+2 — Xj 
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• Similarly, after the {k - 1 )st divided differences, 

f[Xi, X/+1 , X/+2, . . . , Xi+k-^] and f[x/+i , X/_|_2, . . . , X/_|_/(_1 , 
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have been determined. 
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• Similarly, after the {k - 1 )st divided differences, 

f[Xi, X/+1 , X/+2, . . . , Xi+k-^] and f[x/+i , X/+2, • • • , X/+/c_i , Xj+k] 
have been determined, the kth divided difference relative to 

f [X/+1 , X/+2, • • • , - f[Xi, X/+1 , ■ . . , X/+;t-1 ] 
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• Similarly, after the {k - 1 )st divided differences, 

f[Xi, X/+1 , X/+2, . . . , Xi+k-^] and f[x/+i , X/_|_2, . . . , X/_|_/(_1 , 

Xi+k] 

have been determined, the kth divided difference relative to 

f[Xij Xi+^ ) • ■ ■ ) Xi+k-^ ) ■^/+/f] 

f [X/+1 , X/+2, • • • , - f[Xi, X/+1 , ■ . . , X/+;t-1 ] 

• The process ends with the single nth divided difference, 

f[Xi,X2,...,X„]-f[Xo,Xi,...,Xn_i] 



f[Xo,X^,...,Xn] 
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• Returning to the interpolating polynomial, we can now use the 
divided difference notation to write: 

ao = fixo) = f[xo] 
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Using the Divided Difference Notation 



• Returning to the interpolating polynomial, we can now use the 
divided difference notation to write: 

ao = fixo) = f[xo] 
^^^fM^ = /[x„,x,l 

Xi — Xq 

• Hence, the interpolating polynomial is 

P„(x) = f[xQ] + f[xQ, xi](x - xo) + a2(x - Xo)(x - Xi ) 
+ • • • + an{x - xo)(x - Xi ) • • • (x - x„_i ) 
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Pn{x) = f[xo] + f[xo, xi](x - xo) + a2(x - Xo)(x - Xi ) 
+ • • • + a„(x - Xo)(x - xi ) • • • (x - Xn-^ ). 

• As might be expected from the evaluation of Sq and a^ , the 
required constants are 

Bk = f[XQ,X^,X2,...,Xk], 

for each k = 0,^,. . . ,n. 
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Pn{x) = f[xo] + f[xo, xi](x - xo) + a2(x - Xo)(x - Xi ) 
+ • • • + a„(x - Xo)(x - xi ) • • • (x - Xn-^ ). 



• As might be expected from the evaluation of ao and a-\ , the 
required constants are 

Bk = f[XQ,X^,X2,...,Xk], 

for each k = 0,^,. . . ,n. 

• So Pn{x) can be rewritten in a form called Newton's 
Divided-Difference: 

n 

Pn{x) = f[XQ\ + ^ f[Xo, Xi , . . . , X/,](X - Xo) • • • (X - X/f„i ) 
k=^ 
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Questions? 



Reference Material 



Forward Difference Operator A 



For a given sequence {Pn}'^=Q, tlie forward difference Apn (read "delta 
Pn") is defined by 

Apn = Pn+\ - Pn, for /I > 0. 

Higlier powers of tine operator A are defined recursively by 
A''pn = AiA^-^Pn), for/f>2. 



< Return to the Divided Difference Notation 



